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Abstract 

By assuming a self-similar structure for Kelvin waves along vortex loops with successive 
smaller scale features, we model the fractal dimension of a superfluid vortex tangle in the 
zero temperature limit. Our model assumes that at each step the total energy of the 
vortices is conserved, but the total length can change. We obtain a relation between 
the fractal dimension and the exponent describing how the vortex energy per unit length 
changes with the length scale. This relation does not depend on the specific model, and 
shows that if smaller length scales make a decreasing relative contribution to the energy 
per unit length of vortex lines, the fractal dimension will be higher than unity. Finally, for 
the sake of more concrete illustration, we relate the fractal dimension of the tangle to the 
scaling exponents of amplitude and wavelength of a cascade of Kelvin waves. 

1 Introduction 



Turbulence in helium 11, or superfluid turbulence, consists of a tangle of quantized vortex lines 
[H [2]. Until recently, in most experiments superfluid turbulence was created in superfluid 
helium at rest in the presence of a heat flux, the so-called " counterflow" [3l H] , an interesting 
problem of non-equilibrium physics [5^, '6] . More recently, superfluid turbulence was generated 
by agitating the liquid helium using grids or propellers [3 El [9]. Particularly interesting is 
the case in which the temperature T is small enough (T < IK) that the normal fluid fraction 
of helium 11 is negligible, hence viscous dissipation and mutual friction play no role. In this 
low temperature limit, superfluid turbulence takes its purest form: a tangle of reconnecting 
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vortex filaments which move under the velocity field of each other. The importance of vortex 
reconnections, first recognized by Schwarz and later proved by Koplik and Levine [llj . 
cannot be underestimated [12^ I13j . Vortex reconnections randomize the vortex tangle and 
initiate the physical mechanisms of the decay of the tangle's kinetic energy in the absence 
of viscous losses. The first mechanism is the direct conversion of energy into sound in the 
form of rarefaction pulses at reconnecting events, as predicted by the condensate nonlinear 
Schroedinger equation model [14j . The second mechanism is a cascade of Kelvin waves of 
shorter and shorter wavelengths [l5]-[23] triggered by vortex reconnections. This process of 
generation of small scales can proceed without significant kinetic energy losses to spatial scales 
which are small enough that the kinetic energy of the highly curved and cusped fragments of 
the vortices is radiated away as sound [23]~l27] (phonon emission), that is to say, ultimately 
kinetic energy becomes heat. In this regime, one expects the vortex tangle to exhibit fractal 
features, if the mentioned processes act in a self-similar way on several orders of spatial lengths. 
The idea that at very low temperatures the energy can be released by vortex reconnections, 
from smaller and smaller structures, hence it increases the twisting and the winding of the 
superfiuid turbulence, was originally suggested a long time ago by Feynman in his pioneering 
article on the applications of quantum mechanics to liquid helium [28j , before we knew about 
fractals or the Kelvin wave cascade, and was explored in detail by Svistunov |15j . 

Here we propose simple geometrical models of the fractal dimension of superfiuid turbu- 
lence, which represent reconnections and interactions between vortex loops and the subsequent 
formation at the next generation of new vortex loops and Kelvin waves on them. The models 
are too simple to be dynamically realistic, but sufficiently appealing for a qualitative under- 
standing of some physical features infiuencing the fractal dimension. We stress that we are not 
attempting to develop a theory of the Kelvin wave cascade based on actual vortex dynamics, 
but we shall move with simpler considerations. Our motivations are the growing interest in 
superfiuid turbulence at very low temperatures [29l |30] , and previous remarks on the fractal 
nature of superfiuid turbulence. In particular we remind the work of Kivotides et al. [HI], who 
numerically determined a fractal dimension larger than unity (but at finite temperature, not in 
the limit of absolute zero which we consider here) , of Nemirovskii et al. [32] (who considered 
the infiuence of the possible fractal dimension of the tangle on the energy spectrum of the 
turbulent velocity field) and of Jou et al. [33] (who proposed an heuristic form of a Vinen's 
generalized equation for the dynamics of a fractal vortex tangle). 

Our aim is to model the fractal dimension of the tangle under the condition of constant 
energy, but separating the scaling behaviour of vortex length and of vortex energy during the 
transfer of line length to smaller and smaller length scales. The underlying physical idea is 
that the energetic contributions of very close parts of the vortex lines may interfere with each 
other, thus leading to a non additive global result for the total energy of the loop. 

First of all, we derive a general relation between the fractal dimension of the hierarchy of 
self-similar vortex loops and the behaviour of the vortex energy per unit length at different 
length scales. Afterwards, in order to be more concrete and explicit, we propose two simple 
models of hierarchies of self-similar loops, whose behaviour mimics in a simplified way the 
features of a cascade of Kelvin waves. 

Our simple models are partially inspired to the well-known /3-model for classical intermit- 
tent turbulence [S]- [37] and include the infiuence of geometrical and energetic aspects on 
the fractal dimension. We are not aware of applications of the /3-model of classical turbulence 
to quantum turbulence. We think that this model can be useful to grasp some qualitative 
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transfer amongst different length scales. Our approach differs from that of Svistunov [E] in 
that it takes a less detailed, less quantitative form, but it allows a simpler and more intuitive 
view of the complicated process in question. 



2 Fractal dimension and behaviour of the vortex energy per 
unit length at different scales 

Our aim is to look for an expression of the fractal dimension Dp the vortex tangle in 
terms of the microscopic properties previously mentioned, namely, vortex length distribution, 
amplitude and wavelengths of Kelvin waves, and energy density per unit length. To obtain 
the fractal dimension Dp of the vortex tangle, we use the standard definition [351 138] 

log(jV„./jVo) 

Df = - hm - — , , , , (2.1) 

n^oo log(/„//o) 

where A^^ is the number of steps along a curve (or the number of objects of a given size) and 
In the length of a single step (or the size of a given object). 

We assume that the tangle can be statistically described as a self-similar hierarchy of loops, 
whose forms will be discussed in Section 3. The generation with n = 1 corresponds to the level 
of the biggest vortices, which become more abundant and smaller for increasing n. We call 
Nn the number of vortices at the ra-th generation. In the size of each loop, and E'^ the energy 
of each loop. 

Although the specific expression for the fractal dimension depends on the details of the 
model, we express the fractal dimension in terms of the energy per unit length at several 
scales. Note that in our analysis Dp is a property of the ensemble of self-similar loops, not 
of a single loop. In fact, the individual loops are assumed to be regular lines, and not fractal 
lines. Thus, our fractal dimension characterizes the self-similarity properties of the tangle and 
not of individual vortex lines. 

Before proposing an explicit model of hierarchies of vortex loops, we relate the fractal 
dimension defined geometrically in ()2.ip with the variation of the energy per unit length at 
different length scales. 

We assume that E'^ oc 1"' , where a' is a constant scaling exponent; this means that the 
energy per unit length is 

E'Jln^l^'-\ (2.2) 

Therefore, if a' > 1 the contribution to the energy per unit length decreases for lower length 
scales (shorter In); the opposite is true for a' < 1. If a' = 1 then the energy per unit length 
is the same at each length scale. In principle, the exponent a' does not depend on the fractal 
dimension (j2.ip . but it becomes related to Dp if we assume the condition of constant total 
energy at the different loop generations, as mentioned in the introduction. 

According to the previous definitions, the total energy En at the re-th loop generation is 
given by £"„ = NnE'^. Then the condition that the total energy is independent of n can be 
expressed by 

En = En+i hence NJ^ = Nn+ilf^+i- (2.3) 
If n is large enough, equation p.ip implies that a Nn ^^^^ , thus equation p.3p leads to 

^i-(a7i5^)=^^i-K/fl^)_ (2.4) 
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In order that this equahty is true for any n, one must have that Dp = a' . This result shows 
the strong connection between the energetic features of the tangle and its geometrical structure, 
independently of the detailed form of the loops in the hierarchy. When the contribution to 
the energy per unit length of the smaller length scales is smaller than the contribution of the 
larger scales then Dp > whereas in the opposite case Dp < 1. The case Dp < 1 seems 
physically unacceptable because it would imply that vortex line fragments in objects perhaps 
similar to the Cantor's dust, which would violate the condition that the vorticity is solenoidal. 
This result would be supported by the numerical simulation of Kivotides et al. [3l]. The 
only acceptable situation is that the larger length scales contribute more to the energy per 
unit length than the smaller length scales. In view of the meaning of a', we write the fractal 
dimension (j2.ip as 

n^oo log(/nj 

which, by using /„ oc Nn ^^^^ , can also be written as 

Dp n^oo log(iV„) ^ ^ 

In the next section, we shall introduce a model of loop generations which relates the fractal 
dimension to the amplitude and the wavelength of the Kelvin waves. 



3 Geometric and energetic assumptions 

We assume that, as the mutual interaction of vortices induces the formation of helical Kelvin 
waves along the vortex line which also undergo breaking and reconnection processes, the tangle 
can be described as an ensemble of self-similar objects. Neglecting boundaries, we assume that 
these objects are closed vortex loops. Svistunov [15j has also considered the same point of 
view, but with different transformation rules than those we consider here, as we shall discuss 
below. 

The loops can be entangled among themselves in complex topological ways [39j. Here, 
we focus our attention only on geometrical properties such as vortex length, vortex number, 
amplitude and wave-number of Kelvin waves, and energy per unit length. We do not consider 
the topological details of the vortex entanglement. 

We envisage that the generation of vortex loops takes place according to the following 
rules: 

1. We take as reference configuration a collection of A^o vortex loops of length Lq = 27ri?o, 
where Rq is the average curvature radius of the loop, along which an helical structure 
of Nq helical turns lies, all turns being of radius Rq. This structure models in a simple 
way the formation of Kelvin helical waves along vortices, where Rq is the amplitude and 
/iQ = Lq/N'q is the wavelength of waves. 

2. The next generation is assumed to be composed of Nqv smaller vortex loops (where r is a 
multiplication parameter) of length Lq/P (where (3 is another parameter). The following 
generation consists of A'^o^^ loops of lengths Lo//3^, and so on. Thus, the re-th generation 
is composed of Nn = v'^Nq vortex loops of lengths L„ = Lq/(3'^ (see Figured]). Note that 
the generation of smaller vortex loops at the n-th generation comes from the interactions 
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Figure 1: Each vortex loop at n-th generation with average curvature radius i?„ generates 
r = 2 vortex loops at (n + l)-th generation, and so on at the next generation. In the figure a 
constant number of helical turns = Nq (q = 0) is assumed, while the radius R'^ decreases 
according to the rule ii^ = R'q/(3^''^. 

and reconnections among vortex loops of (n — l)-th generation (the dots in each loop of 
Figured] denote the reconnection points). The parameters r and /3 thus depend on the 
details of the dynamics. 

3. Besides the above set of rules which determine the number of loops in the n-th 
generation and the average radius Rn of the main circle of the loop, it is necessary to 
give a second set of rules for and R'^, which are the number of helical turns and 
the average radius of each turn. We interpret R'^ as the amplitude of a Kelvin wave 
and hn = Ln/N'^ as the wavelength. We assume that each loop has N'^ = NQ{r')"' 
helical turns, and that the radius of the helical structure at the n-th generation scales 
as R'^ = R'qI j3"^^ i.e. it scales in a different proportion than the curvature radius of the 
total loop. We introduce two scaling coefficients a and 7, setting r' = and j3' = 
(in Figure [H for the sake of simplicity, it is assumed that the number of helical turns in 
each vortex loop is constant, that is a = 0). 

4. In the breaking and reconnection processes the total energy remains unchanged because 
at the considered length scales there is no friction, hence no energy dissipation. Thus, 
we impose that E^+x = En for all values of n. This assumption has a limiting length 
scale because at some very small length scale sound radiation becomes relevant and the 
fractalization process stops. 

We stress that the first three rules are not supposed to model actual dynamical processes; 
they are only meant to explore the possible consequences of two physical processes which, in 
this context of superfluid turbulence at very low temperatures, are still poorly understood: 
a direct cascade, leading from bigger to smaller loops, and an inverse cascade, leading from 
smaller to bigger loops. Concerning the fourth step, another plausible choice could be the 
invariance of the total length, instead on the invariance of the energy; this provides another 
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way to examine the fractal dimension which leads to Dp = 1. This result follows from the 
general arguments of the Section 2. The assumption that the total length of the vortex line is 
the same at each step n means 

iV„/„ = A^„+i/„+i. (3.1) 
Equation (|2.ip implies that In oc Nn ^^^^ for n high enough, so equation (j3.ip becomes 

(3.2) 

Since relation ()3.2p is true for every n high enough, then Dp = 1. In our opinion it is an 
interesting result because it shows that the fractal dimension is 1 if the whole vortex line 
length is kept constant, that is if vortices themselves do not contribute to lengthen o shorten 
of the total vortex length. But, as pointed out along this paper, in the low temperature limit 
it seems more plausible the invariance of the total energy with respect to the invariance of the 
total length. 

Since the total number of vortex loops at the n-th generation is A^^^, the total energy stored 
on the vortex loops of the ra-th generation will be 

En = NnE'^ (3.3) 

where E'^ is the energy of the loop at the n-th generation. The length of a single vortex in the 
n-th generation is 



L = 2.K«;^1 + (^)^ 2.iVX^^l + (^)' (i?ri)=" ,3.4) 

Finally, we need to specify the energy of each loop. Of course, the energy of a loop depends 
on the helical structure that wraps the unperturbed loop of average radius that is, on 
the number of helical turns N!^ or, equivalently, on the pass of helices. Unfortunately we are 
unable to calculate this energy analytically. The recent work of Maggioni et al. [40J shows that, 
even for the simpler case of a non fractal, single vortex filament, the numerical calculation of 
this energy is difficult, as it converges very slowly. For this reason, we propose three possible 
approximate scenarios for the energy E'^ of a loop at the n-th generation. 

1. The first scenario assumes that E'^ is proportional to the length of the circular axis of 
the helical loop: 

K = P^Ln [log f - l.eisl , (3.5) 

47r L V «o / 

where ps is the mass density of the superfluid component, k is the quantum of vorticity 
(k = h/m, with h Planck's constant and m the mass of the helium atom), and oq is the 
radius of the core of the vortices (of the order of the atomic radius) [Ij . 

2. The second scenario is to assume that for large- amplitude helical turns close to each 
other {R'^ ^> /i„) the helical loop may be considered as a solenoid, and that its internal 
energy E'^ is of the order of its volume, 27ri?„7ri?^ , times the density of the kinetic 
energy of the superfluid induced by the polarized coil. Using for the induced velocity 
an expression analogous to that for the magnetic field in a solenoid we can write for the 
induced velocity Vsi ~ k//i„ and obtain 

K = vrV.^'^n(f^)'. (3.6) 
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3. The third scenario is a more flexible prescription, incorporating also the features of the 
helical structure. The total length of the deformed circle, for n sufficiently high, is 



1 + 



hr. 



(3.7) 



with hn = Lfi/N^ being the pass of the helices at the n-th generation, which can be 
interpreted as the wavelength of the Kelvin wave, whose amplitude is R!^. We assume 
that the whole helical length contributes to the energy, but multiplied by a dimensionless 
factor depending on the length scale, as represented, for instance, by R'^, which modulates 
the relative influence of the helical turns on the energy of the loop. Thus, instead of 
(13.51). we take 



eL 



Ps'^'^ AT' R' f ^" 




logl^ 

ao 



(3.8) 



where 5 is a constant of the order of 1.6. Here the ratio {R'„/Rq)^ with x > ensures 
that the smaller the length scale is (i.e. the smaller R'^ is), the smaller is the contribution 
to the energy; if x < 0, smaller lengths scales have larger contributions to the energy, 
and X = indicates that the energy is proportional to the loop length. Of course, more 
complicated models could be assumed instead of {R'^/Rq)^. The actual value of x should 
be obtained from a first principles calculation of the energy of helical loops of different 
radii and with different separations between successive helical turns, but, as we said 
before, this is very complicated and goes beyond the simple, limited task which we set 
at this early stage of investigation. 

Since Nn/No = r" and In/^o = ^n/2vri?0; when Iq is assumed to be 2ttRq, the fractal dimension 
(12.11) becomes 



Dp 



lim 

n— >oo 



nlogr 



log {N(,R'/Ro) (r""//37") 1 + [Ro/{R'N(,)f (/37-i/r°) 



2n 



1/2 



(3.9) 



Note that here the limit is taken keeping in mind that, when n becomes large enough, the 
radius of the loop Rn cannot be smaller, or of the same order of the vortex core radius oq. We 
distinguish essentially two cases od physical interest: the long wavelength limit {hn ^ Rn) 
and the large amplitude limit (/i„ <C Rn), corresponding respectively to Kelvin waves whose 
wavelengths are larger or smaller than their amplitudes. 



Long wavelength limit ^ R'^ or f3"' ^ > r") 

In the limit of Kelvin waves with small amplitude and large wavelength, the physically most 
plausible scenarios are the first (j3.5p and the third (|3.8p . 

In the first scenario (13. Sp . under the condition -En+i = En, leads to 



[(nln/3 + a)] = - [((n + 1) In /3 + a)] 



(3.10) 



which is true for any value of n only if r = /3 needs. It follows that the condition hn ^ Rn, or 
^7-1 > g^jj read in terms of 7 and q as 7 — a > 1, and the value of the fractal dimension 
can be obtained from (13. 9p : 

Df = 1 if 7-a>l. (3.11) 
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It is reasonable that Dp = 1 because in this scenario the interference between neighboring 
hehcal turns tends to vanish. 

In the third scenario, multiplying (13. Sp times A'^ and requiring that E^^\ — -^n; i.e., 
Nn+iE'^,i = NnE'^, we are lead to 



,(a+l)(n+l) 




[-(n + l)7log/3 + a] 



-n7 log f3 + a] 



(3.12) 



The second term under the square root in p.l2p is dominant and the relation between r and 
13 becomes 

r = (3.13) 

Substituting this relation into I3^~^ /r" > 1 one gets (7 — a — a^x ~ !)/((! + 7X)) > 0- 
Using (fO) and Nn/No = we have 



= 1 + X7 > 1 



if 



x>o. 



(3.14) 



Note that a negative value of x implies Dp < 1, and that the value of x cannot be less than 
— 1/7 > —1. The result that 7 > 1 comes from the relation (7 — a — 07% — 1)/((1 + 7X)) > 0. 
The conclusion that Dp < 1 is not physically reasonable: it would imply that the vortex 
tangle becomes similar to a Cantor's dust and would violate the solenoid condition (a vortex 
is a closed loop or terminates on boundaries, but here we have no boundaries). 

For the sake of completeness, we also consider scenario two ()3.6p although this expression 
of the energy seem to be physically inadmissible in the long wavelength limit. The constraint 
En+i = En leads to r^+^a ^ ^27-1 and 7 - a - 1 > 0. Substitution into ([33]) yields 



Dp 



(as Dp — 1 is positive when 7 > a + 1). 



27 - 1 
2a + 1 



> 1, 



(3.15) 



Large amplitude limit <^ R'^ or < r") 

Suppose that Kelvin waves have amplitude larger than the wavelength. This means that the 
helical turns which wrap loops at the n-th generation must have radius R'^ larger than their 
pass hn- The most plausible assumptions for the energy seems scenario (|3.6p and (|3.8p . 
Using (13. 6p and En+i = En we obtain 



TT 



2n 



r 



TT 



r' AT' ^ ^ 



a \ 2n+2 



(3.16) 



which is valid at any generation n if r^+^" = /3^'^ ^. Substitution into ()3.9p gives 



Dp 



27 - 1 



< 1, 



a + 7 

provided that 7 — a < 1, which, as we have already said, seems physically unplausible. 



(3.17) 
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Using (I3.8P we obtain again (13.12p . but here with the condition 



^7(i+x) =^"+1. (3.18) 

Substitution into the inequality jr"" < 1 yields (7 — a — ajx ~ 1)/(7(1 + x)) < 0- 
Then, substituting (j3.18p into (j3.9p . we obtain 

Df = ^^, (3.19) 
I -ax 



which implies that 



making apparent that 



^ = 4il±^ (3.20) 



Dp>l if x>0, (3.21) 



(otherwise Dp < 1). The last conclusion requires — 1 < x < 0) because < l/Dp < (1 + 
«)/(7(l + x)). 

Again, although scenario one seems to be unphysical, we remark for the sake of complete- 
ness that the condition r = (3, obtained below equation (|3.1Up . leads to 



Dp = > 1 (3.22) 

7 — a 



provided that 7 — a < 1. 



4 Concluding remarks 

Finally, the following two comments are worthwhile. The first is that, at sufficiently low 
temperatures, the energy-conserving process that breaks or lengthen vortices does not continue 
indefinitely, but terminates at sufficiently small scales, where a significant amount of energy 
is dissipated as sound. The dependence of this energy radiation upon the length scale has 
been studied by Vinen [23], and Kozik and B. Svistunov |25j-|27]. According to the Vinen's 
analysis, sound radiation becomes relevant at length scales of the order of Imin — {i^^/^y^^, 
where k is the quantum of circulation and e is the energy communicated to the system per 
unit volume and time, which is proportional to L^. Thus, Imin is proportional to L~^/^. Thus 
sound emission limits the Kelvin wave cascade process considered here. In classical turbulence, 
viscous dissipation plays a similar role, and determines the smallest scale Idiss ^ [z^^/e]^/^ for 
which the celebrated Kolmogorov scaling is valid |34)-[35). 

Second, the correlation between Kelvin waves has not been considered in this simple model; 
in |32j . however, it has been argued it could play a significant role in the fractal properties; 
it would be interesting to explore their contribution in a more detailed model. Our model is 
too simple to do so; it must be stressed that reconnections play a decisive role in the breaking 
bigger loops into smaller loops; without reconnections, energy and momentum conservation 
would forbid this cascade [15] . 

In summary, we have proposed some toy models which allow to interpret the fractal di- 
mension of a vortex tangle in energetic terms. Their energy is not proportional to the vortex 
length — because of the mutual interference of very close parts of the vortex line — and en- 
ergy, rather than the length, is conserved in the breaking and recombination of vortices. For 
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example, very recent work by Maggioni et al. [40] has demonstrated that for complex vortex 
structures such as vortex coils and vortex knots, the energy per unit length is not constant; 
a similar effect may occur on vortex filaments in superfluid turbulence. We have determined 
a relation between the fractal dimension and the influence of the smaller length scales on the 
total energy. If this influence is smaller than that of the bigger length scales, the fractal di- 
mension is higher than 1. This result can be understood in an intuitive way, because energy 
restrictions do not restrict the presence of many small and complicated vortex loops, which 
tend to fill a proportion of space higher than a simple geometrical line. In contrast, if smaller 
scales contribute considerably to the energy, energy restrictions limit the formation of these 
scales and vortex loops will be relatively large and simple. Our results show that when small 
length scales contribute relatively less to the energy than the long scales, the fractal dimension 
Dp is larger than 1. The opposite is not true; one could have a fractal dimension higher 
than 1, but with an essentially linear relation between energy and length. The logarithmic 
dependence of the fractal dimension on the behavior of the energy per unit length at different 
scales may allow to obtain reasonably physical result for Dp without knowing in full detail 
the exact form of the energy contribution of loops at different scales. 

Another interesting result is pointed out in Section 3, and it regards what happens if the 
whole vortex line length is kept invariant at each generation n, instead of the total energy. 
The result is that the fractal dimension has to be 1, that is the assumption that each vortex 
does not contribute to the lengthening or shortening of the other vortices means that vortices 
are not fractals. 
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